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Quantum phase transition is one of the most interesting aspects in quantum many-body systems. Recently, 
geometric quantum discord has been introduced to signature the critical behavior of various quantum 
systems. However, it is well-known that topological quantum phase transition can not be described by the 
conventional Landau's symmetry breaking theory, and thus it is unknown that whether previous study can 
be applicable in this case. Here, we study the topological quantum phase transition in Kitaev's ID p-wave 
spinless quantum wire model in terms of its ground state geometric quantum discord. The derivative of 
geometric quantum discord is nonanalytic at the critical point, in both zero temperature and finite 
temperature cases. The scaling behavior and the universality are verified numerically. Therefore, our results 
clearly show that all the key ingredients of the topological phase transition can be captured by the nearest 
neighbor and long-range geometric quantum discord. 



In recent years topological phases have been intensively studied in condensed-matter systems. Their under- 
standing is relevant to topological quantum computation, which provides the paradigm to store and manip- 
ulate information in topologically nontrivial systems. The topological states are immune to local noise due to 
their nonlocal topological nature 1,2 . The progress was impeded by the fact that non-Abelian anyons are known to 
appear inp-wave superconductors 3,4 . Recently, it has been realized that thep-wave pairing can be emulated with s- 
wave pairing provided with spin-orbit coupling and Zeeman splitting 5,6 . Especially, this provide the possibility of 
direct simulation 7,8 of Kitaev's ID p-wave spinless quantum wire model 4 , which lead to elementary experimental 
evidence 9 for the existence of Majorana fermions. Therefore, these progresses have greatly advanced the field. In 
addition, The analytic eigenspectra in the ID Kitaev model have been obtained, and it was shown that this model 
has one gapless topological phase. But the characters of the topological phase from quantum information 
approach have never been well studied previously. 

Quantum phase transition (QPT), originates from the quantum fluctuations at zero temperature, is character- 
ized by nonanalytical changes in the physical properties of the ground state of a many-body system governed by 
the variation of a parameter of the system's Hamiltonian 10 . Quantum discord (QD) 11 is defined by the difference 
between two quantum analogues of two classically equivalent expressions of mutual information. The QD has 
been introduced to analyze the QPT of the spin XY model 12 and much effort has been devoted to various quantum 
critical systems 13,14 . It is worth to note that other approaches to characterize QPT have also been put forward, such 
as entanglement 15,16 , quantum fidelity 17 , fidelity susceptibility 18 , geometric phase 19,20 and Loschmidt Echo 21 . Since 
the calculation of quantum discord is based on numerical maximization procedure, it does not guarantee exact 
results and in the literature there are few analytical expressions including special cases. To avoid this difficulty, 
Dakic et al. 22 introduced geometric measure of quantum discord (GMQD) which measures the quantum correla- 
tions through the minimum Hilbert-Schmidt distance between the given state and zero discord state. Although 
quantum discord around a QPT for a fermionic lattice hamiltonian has been studied, they only consider the 
magnetic phase transition in the transformed picture instead of the physical picture. The Kitaev model is related 
by the Jordan- Wigner transformation to the transverse field Ising model, however, the Jordan-Wigner trans- 
formation is fundamentally non local, the information is not robust to local perturbations in the experimental 
system. This is because the Jordan-Wigner transformed local perturbations become exactly the non-local sort of 
term that can mess up topologically protected states. The geometric quantum discord of the ground state is 
naturally expected to shed some light on the understanding of topological QPT. 

Traditionally, QPT can be understood in the frame of Landau's continuous phase transitions paradigm with 
local order parameter and long range correlations. However, it has been found that the local order parameters 
cannot describe all possible orders 23 . Therefore, it is interesting to investigate whether this success, describe QPT 
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Figure 1 | GQD D G of the ground state (a) and its derivatives dDJdt (b) 
as a function of t around the critical point f c = 1 for different 
superconducting gap A. The lattice size L = 1001. 

with geometric quantum discord (GQD) D G , can be extended to the 
purely topological QPT cases. Here, we explore this question taking 
the topological QPT in Kitaev's quantum wire model 4 as a example. 
The present work aims to provide a more general understanding of 
geometric quantum discord in topological phase transitions. The 
reason why geometric quantum discord can describe topological 
phases or determine their boundaries is due to the fact that they 
depend only on the properties of the ground state of the system. 
We analyze zero temperature and finite-temperature scaling para- 
meters, extracted from the D G of ground state. First, it is found that, 
in the thermodynamic limit, the derivative of ground-state D G is 
nonanalytic at the critical point. A finite size scaling analysis is car- 
ried out for model with different lattice sites, and the scaling behavior 
is confirmed numerically. This allows us to extract the correlation 
length critical exponent. To confirm the universality, we also per- 
form a finite-temperature scaling analysis for this model, aiming at 
extracting the correlation length critical exponent from the scaling 
behaviors of thermal-state D G . The finite-temperature scaling is 
manifested when the temperature approaches zero, in analogy with 
finite-size scaling. 

Results 

The model and topological phase transition. We consider the 
typical lattice model of the ID p-wave superconductor, which is 
described by the L-site Hamiltonian 



L L-l 

H=^^-^(tctc, +1 +Ac|c] +1 +?i.c.J. (1) 

i — 1 i — 1 

where /( is the on-site chemical potential, t is the nearest-neighbor 
hopping amplitude, A is thep-wave superconducting gap, and c\ and 
Cj are fermionic operators that satisfy the anticommutation relations 
{cj,c m } = <5„ m ,{c„,c m } = {c^,c^} =0. The model has two distinct 
phases, i.e., topologically trivial and nontrivial phases. Under open 
boundary condition, the two phases can be distinguished by the 
presence or absence of zero-energy Majorana bound states at the 
ends. It is well known that the topologically non-trivial phase of 
the model is best illustrated for the choice of parametersfor u/t < 2 
for any Without loss of generality, we assume t and A are both 

real and set ^ = 2. Therefore, our choice makes the model possess a 
critical point at t = 1 for any A?^0. 

GQD close to the phase transition. To illustrate the intrinsic relation 
between the geometric quantum discord of ground-state and 
quantum phase transition in this model. We plot GQD D G and its 
derivative dD G /dt as a function of the Hamiltonian parameters t. As 
shown in Fig. 1(a), given the value of A, the GQD increases with 
increasing the nearest- neighbor hopping strength t. Meanwhile, it is 
shown in Fig. 1(b) that the derivatives of GQD (dD G /dt) exhibit sharp 
peaks around the critical points t c = 1. First, these sharp peaks 
manifest the nonanalytic property of the GQD at the critical point. 
Second, this also verified the assertion made by kitaev 4 that the critical 
point will not change with the changing of A. This means the GQD of 
the ground state can signature the criticality in this model, as in the 
conventional QPT cases. 

Finite size scaling. In quantum critical phenomena, the most impor- 
tant themes are scaling and universality. In a finite lattice model, it is 
expected that the anomalies will become clearer and clearer as the 
size of the lattice increases. The relevant study is the so-called finite- 
size scaling. The critical features are characterized in term of a critical 
exponent v defined by t, = (t — t c )~ v with q representing the 
correlation length. To further understand the relation between 
GQD and quantum criticality, we calculate the derivative of GQD 
dD G /dt for A = 1 and different lattice sizes around the critical point, 
as shown in Fig. 2(a). For simplicity, we first look at the case of A = 1, 
and we will discuss the properties of the family of A 1 later. In fact, 
there is no real divergence for finite L, but the curves exhibit marked 
anomalies and the height of which increases with lattice size. The 
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position of the peak t m can be regarded as a pseudo-critical point 
which changes and tends as L~ 146 towards the critical point and 
clearly approaches t c as L — » o°. This scaling behavior of dD G /dt is 
also clearly shown in in Fig. 2(b). 

We now further going to deal with the critical exponent that 
governs divergence of the correlation length. As shown in Fig. 3(a), 
the singular behavior of dD G ldl for the infinite chain can be analyzed 
in the vicinity of the quantum criticality, and we find the asymptotic 
behavior as 



dDc 
dt 



~Ki ln|r — r c | + const., 



where Ki = —0.0638. On the other hand, As shown in Fig. 3(b), the 
value of dD G /dt at the point t m diverges logarithmically with increas- 
ing lattice size as: 

dDc 
dt 1 



= k 2 In L + const.. 



(3) 



where k 2 = 0.0601. According to the scaling ansatz, the ratio |ki/k 2 | 
gives the exponent v. Therefore, v ~ 1 is obtained in our numerical 
calculation for the lDp-wave superconductor model. Furthermore, 
by proper scaling and taking into account the distance of the max- 
imum of D G from the critical point, it is possible to make all the data 



for the value of F = 10 j^l — exp [dD G /dt — dD G jdt\ ^ J J as a func- 
tion of L 1/l (r — t c ) for different L collapse onto a single curve. The 
result for several typical lattice sizes is shown in Fig. 3(c), where we 
can also extract the critical exponent v = 1, in the vicinity of the 
critical point. 

Universality. As is well known, ID p-wave superconductor with 
nearest neighbor hopping belongs to the same quantum 
universality class for non-zero A, with the same critical exponents 
v = 1 . To confirm the universality principle in this model, we need to 
check the scaling behaviors for different values of the parameter A. 
For instance, from Fig. 3(a) and Fig. 3(b) we get Ki = 0.0716 and ic 2 = 
— 0.0699 for A = 0.8. Moreover, we also verify that, by proper scaling, 
all data for different L but a specific A and a critical point t c = 1 can 
collapse onto a single curve. The data for A = 0.8 are show in 
Fig. 3(c). We can extract the same critical exponent v = 1 from all 
the above results. 

Extended to finite temperatures. We would like to further study the 
relation between the thermal-state D G and quantum phase 
transitions at a finite temperature. The derivative of GQD dD G /dt 
as a function of t at different temperatures T (including zero 
temperature) are presented in Fig. 4(a). At zero temperature the 




L(t-tm) 

Figure 3 | (a) Plot of dD G /dt against In It — t c \ for evaluating the thermodynamic approaching to the critical point t c = 1. The line slope close to the 
critical point is — 0.0638( — 0.0716) for A = 1.0 (A = 0.8). (b) The maximum value of the derivative dDc,/dt\ t=t at the pseudocritical point t m as a 
function of lattice sizes. The line slope close to the critical point t c = 1 is 0.0601 (0.0699) for A = 1.0 (A = 0.8). The critical exponent v for 
the correlation length is determined by the two slopes in (a) and (b) for a fixed A. (c) A finite size scaling analysis is carried out for the value of 
F= 10 3 ^1 — exTp^dDc/dt — dD^j dt\ t=t \ as a function of L(t — t m ) for different lattice sizes. All the data from L = 51 up to L = 401 for a fixed 
parameter A collapse on a single curve. This shows that the system at the critical point is scaling invariant and that the correlation length critical exponent 
v = 1. 
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derivative of D G shows a singularity at t c = 1, but at nonzero 
temperature, there are no real divergence. Nevertheless, there are 
clear anomalies at low temperature, and the height of which 
increases with the decrease of the temperature. This can be 
regarded as the precursors of the QPT. What is more, the position 
of the maximum derivative (pseudocritical point) t m changes and 
tends as T 144 and clearly approaches t c when T — » 0, as shown in 
Fig. 4(b). Meanwhile, the maximum value of dD G /dt at the 
pseudocritical point t m diverges logarithmically with the decrease 
of the temperature 

dDg 

~-r±\ tm *=>K 3 ]nT+ const. (4) 

Our numerical results, as shown in Fig. 4(c), give jc 3 = 0.0611. The 
ratio of the two slopes (k 1 /k 3 ) for a fixed parameter A is equal to the 
critical exponent v = 1. Moreover, we also verify that by proper 
scaling, all data for different temperatures T but a specific A will 
collapse onto the same curve. The data for A = 0.8 and A = 1.0 
near the critical points t c = 1 are shown in Fig. 4(d). The numerical 
results agree with the finite-size scaling ansatz and the universality of 
the above results. 

It has also been noted that unlike pairwise entanglement, which is 
typically short ranged, quantum discord does not vanish even for 




In T L(t-tm) 

Figure 4 | (a) The derivatives of GQD dDJdt as a function of the 
parameter t. The curve corresponds to different temperatures T = 0, 0.01, 
0.02, and 0.04. With the decrease of the temperature, the maximum gets 
pronounced, and the pseudopoint f,„ changes and tends as V M towards the 
critical point t c = 1 (b). (c)The maximum value of dDJdtzX. the 
pseudocritical point f,„ of this model as a function of temperature T. The 
slope of the line is -0.0611 (-0.0717) for A = 1(A = 0.8). (d) The value of 
F=10 2 (l-exp(dD G /dt-dD G /dt\ t=t J) as a function of ( t - tJ/Tfor 
different temperature. 



distant lattice pairs. Fig. 5 shows that quantum discord provides the 
expected long-range behavior of quantum correlations LD G , L = 15 
for ID Kitaev model exhibiting QPTs. In a topological nontrivial 
case, one can observe stronger quantum discord. The critical expo- 
nent v = 1 for the correlation length is determined by the two slopes 
in Fig. 5(c) and Fig. 5(d) for a fixed A. Such a behavior close to 
topological QPT is remarkable, even though quantum discord has 
been shown to be non-vanishing. Concerning QPTs, long-range 
quantum discord has been considered as an indicator of the critical 
point of topological QPT in several systems, having succeeded in this 
task even in situations where entanglement fails, which may shed 
some light on the understanding of topological phase transition for 
the other more complex systems from the view of quantum discord. 
Finally further numerical calculation of the quantum discord shows 
that for 2D Kitaev model exhibiting topological QPTs in the range of 
t x + t y = 1, critical exponent v also gives exactly the same results with 
the theoretical results z = 2, v = 1/2. 

Discussion 

Topological states that are protected from local perturbations need 
to be supported by quantitative calculations, in particular working 
with quantum systems of finite sizes. To study the experimental 
feasibility, we add a random defect perturbation to each hopping 
amplitude to check the robustness of the topological QPT under the 
perturbation. Taking the system defects with parameters as given in 
Fig. 6, the numerical results indicated that the scaling and the 
critical point of topological QPT are robust against defects. We 
find topological QPT located in the boundary is robust even the 
defect amplitude takes 8 = 0.1. In summary, we have performed a 
finite size scaling analysis, whose analytical expression has been 
extracted from the D G corresponding to different values of the 
control parameter near the critical point. This makes it possible 
to extract the correlation length critical exponent. Finally, our study 
establishes the connection between the D G and topological QPT at 
both zero temperature and nonzero temperature. All key features of 
the quantum criticality, such as scaling, critical exponent, the uni- 
versality are presented. We also would like to point out that the 
results obtained in this paper does not depend on the model. An 
interesting question is that when added with next-neighbor tunnel- 
ing, i.e., H 2 = X 2 ^cJc,+2 + Ac,c 2 +h.cA , new critical points will 
appear and this needs further study. 

Added note. In the revise manuscript we find for the same question, 
Luo etal 24 also investigate the fidelity susceptibility (FS) and the 
topological phase transition in a two-dimensional spin-orbit 
coupled (SOC) Fermi superfluid. 

Methods 

Diagonallzation of the model hamiltonian. Under the open boundary condition, 
the Hamiltonian can be immediately diagonalized 

I 

H=I>0)jA-l/2), (5) 

H = l 

by introducing a Bogoliubov-de Gennes (BdG) transformation 

fj n — ^^ _ 0 (wn.iCj +Vn i( cjV where ). n denotes the quasi-particle energy, and u nA and 

v nii can be obtained by solving the corresponding BdG equations, the Bogoliubovde 

Gennes matrix is written in blocks as Hbdg — I - } with 

V-A -hj 

hjj = 8ijH i — t8ij+j — and = A — c^+iA. In order to obtain the D G in 

the following, we need to calculate the ground-state correlation functions 

C?j° (cjCj), C-/ - hicj\ C)f = (clcj\ and C)j = ($$\ written in terms of the 
amplitude of the BdG transformation as 
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Figure 5 | (a)(b) Long-range quantum discord LD G and its derivatives dLDJdt as a function of the parameter f. (c) dLDJdt against In ] t — t c \ 
for evaluating the thermodynamic approaching to the critical point t c = 1. (d) The maximum value of the derivative dLDc/dt\ t=t „ at the 
pseudocritical point t,„ as a function of lattice sizes. The critical exponent v for the correlation length is determined by the two slopes in (c) and (d) for a 
fixed A. 
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»=1 n = l 



Computation of D G . An arbitrary two-qubit state can be written in Bloch 
representation: 



(6) 



1+^2 (xiCi<2>i+yii®Oi) + ^2 RijGi®Cj 



(7) 



where — Trp(<7j ® /), y { = Trp(7 (x) a-;) are components of the local Bloch vectors, 
(Tp i e {1, 2, 3} are the three Pauli matrices, and — Trp(Oi ® Oy) are components of 
the correlation tensor. For two-qubit case, Then a analytic expression of the D G is 
given by 22 



D G (p) = i (|W| 2 + P|| 2 -M 



(8) 



where x = (x u x 2 , x 3 ) T and fc max is the largest eigenvalue of matrix K = xx T + RR T . 
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